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ABSTRACT 

Understanding the cosmic re-ionization is one of the key goals of the modern observational 
cosmology. High redshift QSO spectra can be used as background light sources to measure 
absorption by intervening neutral hydrogen. We investigate neutral hydrogen absorption in a 
deep, moderate-resolution Keck/Deimos spectrum of QSO CFHQSJ2329-0301 at z=6.4. This 
QSO is one of the highest redshift QSOs presently known at z=6.4 but is 2.5 mag fainter 
than a previously well-studied QSO SDSSJl 148H-5251 at z=6.4. Therefore, it has a smaller 
Stromgren sphere, and allows us to probe the highest redshift hydrogen absorption to date. 
The average transmitted flux at 5.915 < Zabs < 6.365 (200 comoving Mpc) is consistent with 
zero, in Ly a, Ly /?, and Ly 7 absorption measurements. This corresponds to the lower limit 
of optical depth, rfy^ > 4.9. These results are consistent with strong evolution of the optical 
depth at z > 5.7. 

Key words: quasars: individual, cosmology :early universe, black hole physics, galaxies: high- 
redshift 
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1 INTRODUCTION 

Understanding when and how the neutral dark Universe was re- 
ionized by the first stars, galaxies, and black holes is a key goal 
of modem observational cosmology. The bright UV continuum of 
high-redshift QSOs can be used as a background light source to 
probe this cosmic re -inonizatio n since it is absorbed by neutral hy- 
drogen. Previously, iFan et al] (2006, F06 hereafter) compiled 19 
QSOs at 5.74 ^ 2; < 6.42 to find the sudden increase of optical 
depth (r) and the length of dark gaps at z~6, suggesting an ap- 
proach to the end of the re-ionization. 

However, increasing diversity has been observed in the optical 
depth in various line of sights, suggesting it is important to use a 
statistically significant sample size to discuss the evolution of opti- 
cal depth. The highest redshift QSO previously used for absorption 
analysis, S PSS J 114 8-1-5251 at z=6.42, had escaping flux detected 
at z=6.25 jpan et aL h2003). This means there is still room for the 
QSO absorption measurement to be a useful probe of the cosmic 
re-ionization. 

In this work, we observe QSO CFHQSJ2329-0301 at z=6.4 
with the Keck/Deimos spectrograph. This QSO is ~2.5 mag 
fainter than the previously well-studied highest redshift QSO 
SDSSJl 148-1-5251 at z=6.4. Therefore, the QSO has a few times 
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smaller Stromgren sphere (in radius), allowing us to measure opti- 
cal depth much closer to the QSO emission at higher redshift. We 
measure the neutral hydrogen absorption in this QSO spectra to 
gain a better understanding of the evolution of the optical depth at 
high redshift. 

Unless otherwise stated, we adopt the WMA P cosmology: 
(h, Qm, ^l) = (0.7, 0.3, 0.7) jKomatsu etalj201ll) . 



2 OBSERVATION 



Our target is a QSO CFHOSJ2329-0301l! willott et al.ll2007 l W07 
hereafter). Initially the redshift was reported to be z=6.43 based on 
Ly a (W07). The redshift was later corrected to be z=6.417 based 
on the Mg II line iWillott et alj201(ih because blueward of the Ly a 
emission line suffers from heavy absorption from neutral hydrogen. 
This is one of the three highest redshift QSOs to date, along with 
SDSSJl 148-1-5251 at z=6.41 9 (Fan et al. 2003) and CFHQS0210- 
0456 at z=6.438 l lWillott et al.i2010l) . This QSO is known to be in a 
dense environment surrounded by 7 Lyman break galaxies (LBGs) 
jUtsumi et al.ll201(]l) and to have extended Ly a emission and con- 
( lGotoetal.i2009l) . Characteristics of the target are summa- 
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rized in Table[T] 

We obtain ed moderate reso lution spectra of the QSO using 
Keck/Deimos tFaber et alj|2003l) on Sep. 12-13, 2010. Although 
Deimos was used in the multi-slit mode to investigate the environ- 
ment of the QSO, we focus on the QSO spectra in this paper We 
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Figure 1. Two dimensional spectra of CFHQSJ2329-0301, taJcen witli tlie Keck/Deimos spectrograpli. One-dimensional spectram is shown in Fig.|2] 




Figure 2. One-dimensional spectra of CFHQSJ2329-0301. The blue dashed Hne shows the best-fit continuum level, assuming a power law of (x v~ . 
The cyan lines show the wavelengths masked out due to the strong sky emissions. The green Unes show the sky spectram in arbitrary units. The gray lines 
show the first three wavelength regions where the flux was consistent with zero. The corresponding two-dimensional spectram is shown in Fig.[T] 



used the 830 lines mm^ grating and the OG550 filter with a cen- 
tral wavelength of 8OOOA . This setting has a wavelength coverage 
of 6OOOA to lOOOOA with a spectral resolution of ~2.5A for 0.75" 
slit. The spatial resolution is 0.1 185" pixen^. The total integration 
time was 5.5 hours. 

W e used the DEEP2 pipeline l lWirth et al.ll2004l : lOavis et al.l 
to reduce the data, except the background subtraction, which 
we did manually to remove ghost features of the 830 grating. Wave- 
length calibration is based on the HeNeAr lamp. The spectrophoto- 
metric standard G191-B2B was observed and used for flux calibra- 
tion. Fig.[T]shows two-dimensional spectrum of the QSO. 

At >7500A, many strong sky OH emission lines are present. 
Although the DEEP2 pipeline subtract these sky lines well, the 
noise level at the exact wavelengths of the sky lines are significantly 
higher than the other wavelength clear of emission. Therefore, we 
have masked out these wavelengths of strong sky emissions from 
the absorption measurement. In Fig.|2] the wavelengths masked out 
are shown in thin cyan lines. The green lines at the bottom show 
the sky spectrum in arbitrary units. We checked that this procedure 
only affects our measurements at the few per cent level, and does 
not change the main conclusion of the paper. 

We used double Gaussians to fit Ly a and NV lines. The 
double-gaussian model do not fit the observed spectra well due to 
the heavy absoiption blueward of Lya and broad natures of both 
lines. The best-fit rest-frame equivalent widths are 110±60A, and 



6±lA for Ly a and NV, respectively. These values are consistent 
with previously reported values (W07). 

The spectrum shows four strong Ly a absoiptions on the blue- 
side of the emission (Fig.|2l(. There also exist absorption lines at 
9183, 9212, and 9240A with equivalent widths of 1.4, 1.5 and 
2.6A, respectively. Because none of the pairs corresponds to dif- 
ferent lines of the same system (such as SiIIA1260, SiIIA1304, 
OIA1302, and CIIA1334), the identities of these lines are not clear. 
There also is an emission line at 9681 A, corresponding to OI line 
(1302.17A)atz=6.43. 



3 ABSORPTION MEASUREMENTS AND RESULTS 

High redshift QSOs can be used as a background light source to 
investigate the inter-galactic medium (IGM) at high redshift. In the 
following subsections, we measure average absorption in the Lya, 
Ly P , and Ly 7 forest regions. 

To measure the IGM absorption, we need to determine the 
continuum before the absorption. This is not a trivial process be- 
cause the continuum blueward of the Lya is heavily absorbed by 
the neutral hydrogen. At the same time, our optical CCD loses sen- 
sitivity at around 10,000 A, and we do not have a long baseline to fit 
continuum, let alone the presence o f strong Lya, NV and 1 emis- 
sions. Following previous work bv iFan et al.l (l2000l l200lh . we as- 
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Table 1. Target information adopted from lWillott et al.l l l2007l) and lGoto et alH2009l) . 



Object 


ZMgll 




J 


CFHQS J232908.28-030158.8 


6.417±0.002 


25.54±0.02 21.165±0.003 21.683±0.007 


21.56±0.25 




Figure 3. Lyman transmission in the spectr um of CFHQSJ2 329-0301 (red, 
filled-circles). T he black t r iangle s are from lFan et al.l ( |2006|) . and the gray 
squares are from lSongaiial | |2004 . 

sume a power-law continuum of oc i^^"'^. When the continuum 
slope is uncertain, the slope of i^""'^ is commonly used by previous 
works. This is the case for most QSOs at z~6. We fit the continuum 
of r/""'^ to the spect rum at > 9300A to det ermine the continuum 
level. We also follow iFan et"^ ( I200(3ll200lh to consider variation 
from z/"*""'^ to u^^'^ to estimate errors of r. However, this is a con- 
servative approach since at least at low er z, most QSO continuum 
slopes are between -0.3 and -0.9 (e.g., ^ Vanden Berk et al" I l200lh . 
Some recent papers do not even include variation of the slope in the 
errors of r. In this work, errors in the continuum fit and the variation 
in the slope are included in the final uncertainty in the transmission. 
Based on this continuum, we measure average transmission in the 
following narrow wavelength (redshift) windows. Then, the trans- 
mission is used to calculate the effective optical depth. 

3.1 Ly a absorption 

With QSO CFHQSJ2329-0301, we can measure absorption closer 
to the Lya emission than bright QSOs, because this QSO has a 
small Stromgren sphere (ionized zone) of 6.4 Mpc (where Lya 
flux drops first to T < 0.1 for spectra binned in 20A pixeP^) 
due to its faint intrinsic luminosity (W07). To measure absorp- 
tion, we used the maximum redshift, Zmax, of 6.365, where blue- 
ward of continuum does not show any sign of rise. This gives us 
Zem — Zmax ~ 0.052. This is the highest redshift absorption mea- 
surement from the QSO spectra to date due to the high redshift of 
QSO (z=6.4) and its small Stromgren sphere. 

We measure Lya absorption down to 1040 A, which is the 
shortest wavelength that is not affected by the Ly/3-l-OVI emissions. 



Figure 4. Effective Gunn-Peterson optical depth from the spectrum 
of CFHQS J232 9-0301 (red, filled-circles). Th e black t r iangle s are from 
lFanetaljj2003) . and the gray squares are from lSongaiial ( |2004|) . The solid 
line is the best power-law fit to the data at 2 < 5.5 by F06 (eq.5). 



Within this range, we used the redshift interval of Az of 0.15, 
which corresponds to ^^60 Mpc in comoving distance, for compar- 
ison with previous work. Although the resolution of our spectrum 
is lower than high-resolution spectra used in the past, this is not a 
disadvantage for a binned measurement like the transmission. The 
first three bins are shown in gray horizontal lines in Fig. |2] The 
median redshifts of five bins are presented in Table |2] Using this 
redshift range and interval, we measure the transmission as 

T{z,ts) =< /""v/r' > 

The measured transmitted flux (T) is shown in Table [2] and 
plotted in Fig. [3] We estimate errors by changing the continuum 
slope from v'^'^''' to v~^'^, quadruply added by the noise in spectra 
and the errors in determining the continuum level. The errors on T 
are dominated by the continuum determination. 

In Fig.[3]and Tabled the first three bins at 5.915 < z < 6.365 
have no detected flux and are consistent with zero within the error. 
Only two lower redshift bins at 5.615 < z < 5.915 have detected 
flux. 

It is conventional to express optical depth r in terms of r = 
—ln{T). We present r as a function of redshift in Fig.|4] For bins 
where we do not have detected flux, we convert the errors on T to 
the lower limit of the r. 

In Figs[3]and|4l sm all triangles and g ray squares represent pre- 
vious data from F06 and lSongailal ( |2004|) . The solid line shows the 
best power-law fit to the data at 2: < 5.5 by F06 (eq.5). 
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Table 2. IGM absorption towards the QSO CFHQSJ2309-0301 at ^ = 
6.42. The widths of bins are Az = 0.15, 0.15, and 0.06 for Ly o , Ly /3 , 
and Ly 7 absoiptions, respectively. 



Redshift 


Line 


Transmission 




6.29 


Ly a 


-0.007±0.008 


>4.9 


6.14 


Ly a 


-0.006±0.006 


>5.1 


5.99 


Ly a 


0.002±0.007 


>4.9 


5.84 


Ly a 


0.073±0.013 


2.6±0.1 


5.69 


Ly a 


0.068±0.013 


2.7±0.2 


6.29 


Ly/3 


0.021±0.041 


>7.2 


6.14 


Ly/3 


-0.002±0.043 


>7.1 


6.29 


Ly7 


0.26±0.23 


>6.4 


6.23 


Ly7 


-0.08±0.23 


>6.5 



3.2 Ly /3 absorption 

At the same neutral hydrogen density, the optical depth r is pro- 
portional to /A, where / and A are the oscillator strength and 
rest-frame wavelength of the transition. Therefore, the r of Ly/? 
is factor of 6.2 smaller than that of Lya. Thus, Ly/3 can probe into 
much more neutral hydrogen than Lya absorption. In this section 
we measure optical depth using Ly/3 absorption from the spectra. 

We assume the same continuum as was used for Lya. We 
choose the minimum redshift of Ly/3 as 970A, above which it is 
not affected by Ly7 absorption. This leaves us the redshift range 
of 6.065 < Zabs < 6.365. The Ly/3 absorptions overlap with Lya 
absorption from lower redshift. Therefore, the Ly a absoiption has 
to be corrected in measuring r^. We used eq. (5) of F06 to esti- 
mate Lya absorption from lower redshift, and corrected the Ly/3 
transmission measurements. 

Table. [2]shows the results. They are also plotted in Fig.|3]with 
the blue diamonds. No flux is detected in Ly/3 measurements. The 
two measured Ly /3 transmissions at 6.065 < Zabs < 6.365 are 
consistent with zero flux within the error. 

To convert Ly/3 transmission to r/3, one has to consider dif- 
ferent optical depths between Ly a and Ly /3 . The difference de- 
pends on the UV background, dumpiness of the IGM, its equation 
of state, and the uniformity of the UV b ackground. According to 
simulations (e.g., lOh & Furlanettoll2005h and empirical measure- 
ments (F06), the Ta/rp conversion is in the range of 2.5-2.9. The 
Ta/r-y conversion lies in the range of 4.4-5.7. Following the discus- 
sion in F06, we use Ta/Tp=2.25 and ra/r^=4.4. 

Fig.|4]shows constraints on the from Ly j3 absorption mea- 
surements. Note that r is converted to Ly a optical depth in Fig.|4] 
Since flux are not detected in both bins. Table |2] shows the lower 
limit in from the error. 

3.3 Ly 7 absorption 

Similarly to Ly /3 , the Ly 7 optical depth is a factor of 17.9 smaller 
than that of Ly a , providing us with a chance to probe into more 
neutral hydrogen. However, the Ly 7 absorption measurement is re- 
stricted by the presence of Ly 5 absorption at the low-redshift end 
and the QSO flux at the high-redshift end. Therefore, we use a 
smaller binsize of Az of 0.06. The Ly 7 transmission measurements 
are corrected for overlapping lower redshift Ly a and Ly /3 absorp- 
tions using the best-fit power laws to lower redshift data (eqs. 5 
and 6 in F06). The optical depth is converted to the r^, using 
rQ/r.^=4.4. The results are in Table|2]and Fig.|4] 



The Ly 7 transmission has more potential since it has a factor 
of 17.9 smaller optical depth. However, the Ly 7 absoiption mea- 
surements are much more challenging. Since the overlapping fore- 
ground Ly a and Ly /3 lines absorb ~98% of the continuum flux in 
the wavelength ranges of Ly 7 , so we need to measure absorption 
in the remaining ~2%. Errors in continuum determination are also 
larger since we use the bluer part of the spectrum. By just changing 
the slope of to u"^'^, we get 10% error in the measurement. 
Combined with the noise in the spectra, we did not obtain tighter 
constraints than Ly /3 . The purple crosses in Fig.|4]show the lower 
limit from the measurement error. 



4 DISCUSSION 

It has often been suggested that QSOs are an ineffective re- 
ionization probe because Lya absorption saturate at neutral hydro- 
gen fractions of /hi ~ 10~*. However the Gunn-Peterson op- 
tical depths of Ly/3 and Ly7 are factors of 6.2 and 17.9 smaller 
than that of Lya. It has also been suggested that metal lines 
such as 01 can be used to probe into even larger optical depths 
( iBecker et alj2006l , l2oT ij). While the dar k gap statistic method can 
also probe into more neutral hydrogen dSongaila & Cowiell2002l : 
IPentericci et al.ll2002l) . our medium resolution spectrum did not al- 
low us to measure the dark gaps. Although we did not detect any 
flux at 5.915 < Zabs < 6.365, another z=6.4 QSO had a flux de- 
tected at z=6.29 with Ly 7 (F06). Therefore, there is room for high- 
z QSOs to be an important re-ionization probe beyond z=6.4. Also 
in both Figs.[3]and|4l r and T show a significant variation, which 
suggests the importance of finding more bright QSOs at z > 6.4 as 
an end of re-ionization probe. On-going and future imaging surveys 
at > 1/im such as UKIDSS, VISTA, CFHQSIR, and PanSTARRS 
are expected to find such QSOs in the near future. 

One of the remaining problems is an uncertainty in continuum 
determination. We used a slope of i/^"'^, but 1^+"'^ and u^^'^ also 
fit continuum at >9300A well. Variation from these is larger than 
the eiTors on the flux, even with our relatively short exposure of 
a faint QSO. Ironically, it is hard to measure intrinsic continuum 
blueward of Lya emission because of the heavy absorption that 
we aim to measure. To make the matter worse, observations of a 
larger sample of QSOs at lower redshift show a continuum break at 
1200< A < 1 300 A, with a stee per power-laws of at shorter 

wavelengths jTelferetal.ll2002h . If there is an evolution, the con- 
tinuum blueward of high-z QSOs remain uncertain. Therefore, the 
continuum determination is a problem in measuring r even if one 
takes very high S/N spectra. 

Interpretation of the evolution of r is a matter of active debate. 
At z > 5.7, F06 reported a departure of r from the best-fit power 
law from z < 5.5, arguing that z ^6 is the end of the reioniza- 
tion, where the Universe was not completely ionized. The average 
length of dark absorption gap also jumps from < 10 to >80 comov- 
ing Mpc at 2 > 5.7. Our r measurements are clearly above the 
solid line (fit to the z < 5.5 data) in Fig.|4] and thus are consis- 
tent with this sudden increase of t at z > 5.7. On the other hand, 
iBecker. Ranch, & Sargenll ( l2007h showed that the probability dis- 
tribution of Ly a transmitted flux of QSOs at 2 < z < 5 are better 
fit by lognormal optical depth distributions. This lognormal fit is 
consistent with z > 5.8 data for no evolution in the hydrogen neu- 
tral fraction. However, W07 pointed out that the lognormal fit does 
not reproduce the data at 5.4 < z < 5.7. 

In either case, the large variation in the r is the primary reason 
why determining a more realistic model is difficult. Although our 
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results add r measurements at the highest redshift probed by QSO 
absorption, it is clear that we need a statistically significant sample 
(>20) of bright QSOs in every A2:=0. 1 at z > 6, where r quickly 
increases, to overcome cosmic variation of t and more accurately 
trace the evolution of r. 



5 SUMMARY 

Using a deep, medium-resolution Keck/Deimos spectra of a QSO at 
z=6.4, we measured neutral hydrogen absorption, and thereby, the 
optical depth. Because our QSO is by 2.5 mag fainer than the previ- 
ously studied QSO SDSSJl 148-1-5251 at z=6.4, we can investigate 
absorption much closer to the QSO emission, reaching z=6.365, the 
highest absorption redshift investigated to date. Our measurement 
of transmitted flux is consistent with zero at 5.915 < 2:"^^ < 6.365 
in all Ly Q , Ly /3 , and Ly 7 measurements. We obtained a lower 
limit of > 4.8. These results are consistent with strong evolu- 
tion of the optical depth at z > 5.7. 
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